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Abstract 

The spectral function related to the correlator of two colour-electric fields along a Polyakov 
loop determines the momentum diffusion coefficient of a heavy quark near rest with respect to 
a heat bath. We compute this spectral function at next-to-leading order, O(a^), in the weak- 
coupling expansion. The high-frequency part of our result (u ^ T), which is shown to be 
temperature-independent, is accurately determined thanks to asymptotic freedom; the low- 
frequency part of our result (w ^ T), in which Hard Thermal Loop resummation is needed 
in order to cure infrared divergences, agrees with a previously determined expression. Our 
result may help to calibrate the overall normalization of a lattice-extracted spectral function 
in a perturbative frequency domain T <^ u <^ 1/a, paving the way for a non-perturbative 
estimate of the momentum diffusion coefficient at w — )• 0. We also evaluate the colour-electric 
Euclidean correlator, which could be directly compared with lattice simulations. As an aside 
we determine the Euclidean correlator in the lattice strong-coupling expansion, showing that 
through a limiting procedure it can in principle be defined also in the confined phase of pure 
Yang-Mills theory, even if a practical measurement could be very noisy there. 
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1. Introduction 



In order to theoretically model various observations made at relativistic heavy ion collision 
experiments, it would be important to quantitatively determine a number of "transport coef- 
ficients" of Quantum Chromodynamics (QCD) at temperatures around a few hundred MeV. 
Such transport coefficients include for instance the shear and bulk viscosities and the thermal 
conductivity, which play a role in the hydrodynamic equations describing the expansion of the 
thermal fireball; as well as heavy flavour diffusion coefficients, the jet quenching parameter, 
and the electrical conductivity, which dictate how various "probes" propagate with respect 
to the hydrodynamically expanding thermal medium (for recent reviews, see refs. [I112])- 

Among all the transport coefficients, the theoretically most accessible one appears to be 
the one related to the diffusive motion of heavy quarks (charm or bottom) near rest with 
respect to the thermal medium. In a classical Langevin picture, the heavy quarks satisfy the 
equation of motion dpi/dt = —rjiypi + ^j, where pi is their momentum and is a random 
force; the transport coefficient % can be referred to as the "kinetic thermalization rate" 
or the "drag coefficient". The simplicity of its determination is due to the fact that it can 
be fluctuation-dissipation related to another coefficient, k = | X^i('?i(^)Ci(0)); called the 
"momentum diffusion coefficient" (the relation reads r/n = K/(2TMkin)(I + C'(T/Mkin)), where 
Mijin is a certain heavy quark mass definition); can in turn be identified as the Lorentz 
force, which for small velocities is proportional to the electric field strength. Thereby k can 
be determined from the two-point correlator of colour-electric fields along a timelike Wilson 
line [31 H] . The correlator of colour-electric fields is simpler to handle theoretically than that 
of the usual bilinear quark currents, because fewer physical scales enter the problem (heavy 
quark related scales like M^in do not appear) and because no transport peak is present. 

Within this framework, a number of results have recently been obtained. First of all, 
the leading perturbative expression for k [5]~IZ] has been supplemented by a next-to-leading 
order correction [8], which was furthermore shown to be so large as to question the valid- 
ity of the weak-coupling expansion. (This computation can be contrasted with the corre- 
sponding results for theories with gravity duals, in particular for strongly-coupled = 4 
Yang-Mills theory in the large- A'^c limit [Si [3], in which case k appears to be much "larger" 
than according to leading-order QCD.) Second, the frequency dependence of the correlator 
k{uj) = 5 e**^* ^j(^j(t)^i(0)) has been determined within various frameworks, such as 
strongly coupled J\f = 4: Yang-Mills theory in the large- limit [31 [lO] , Hard Thermal Loop 
resummed perturbation theory [4], as well as classical lattice gauge theory jTT]. Gaining un- 
derstanding on the frequency dependence is important because a non-perturbative extraction 
of K from lattice QCD would effectively necessitate an ansatz for the shape of k{uj). 

Apart from these theoretical developments, heavy quark diffusion has also been embedded 
in phenomenological models (see, e.g., refs. [71 [12] and references therein). Comparing with 
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experimental data [13J, it appears that the physical value of k would need to be significantly 
larger than the leading-order [5]-[7] or even the next-to- leading order one |8j, perhaps in 
qualitative accord with the AdS/CFT result [9l[3] (cf. also ref. [14J and references therein). 

The purpose of this paper is to address the frequency dependence of ^(a;) within QCD 
and pure Yang-Mills theory, with an eventual lattice determination of k = k(0) in mind. 
To this end, we determine k{ijj) at next-to-leading order, ©(a^), in the weak-coupling ex- 
pansion. Thanks to asymptotic freedom, our result becomes increasingly accurate in the 
large-frequency limit where, as we show, the result is also independent of the temperature. 
In fact, in several respects, the result at (x> > T is analogous to that obtained a long time ago 
for the spectral function related to the electromagnetic current [15]; nevertheless, we are not 
aware of a previous computation of a spectral function at 0{a1) consistently applicable in 

1/2 

the whole frequency range from uj < as T up to cj ^ T. 

The paper is organized as follows. In section [2] we define the basic observable considered and 
introduce our notation. The main steps of the computation are outlined and to some extent 
also carried out in section [Sj a number of technical details are relegated to two appendices. 
Section [4] collects together various results, both in analytic form as well as numerically. 
Section [5] is complementary to the main body of this work, addressing the colour-electric field 
correlator in the lattice strong-coupling expansion. We summarize and offer an outlook in 
section [6l 



2. Basic definitions 

The computation is carried out within continuum QCD, with Nc colours and A'f massless 
quark flavours. The theory is regularized dimensionally, and renormalized quantities are 
expressed in the MS scheme. The only parameter requiring renormalization is the gauge 
coupling; the bare gauge coupling is denoted by gs, the renormalized one by g, and Os = 
g'^/Air. The relation between g^ and g will only be needed at 1-loop order, 

5b = 5 +-^^ 3^ + 0(5), (2.1) 

where is a scale parameter related to "minimal subtraction", and the corresponding MS 
scheme quantity is denoted by p? = 4Tre~'^^ fi'^ . The factor /U~^^ is normally not displayed 
explicitly; the same holds for corrections of 0{e) which are not multiplied by divergent coef- 
ficients. A group theory factor appearing frequently is defined by Cp = {N^ - l)/2Nc. 

The observable considered contains colour-electric fields and temporal Wilson lines. Em- 
ploying sign conventions where the covariant derivative reads = 9^1 — igB^fi, with yl^ 
a traceless and hermitean gauge field, we define a temporal Wilson line at spatial position 
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(2.2) 



r = through 

U{n,Ta) = l + ig^ rdTAo{T,0) + {igj,)^ j 'dr / dr' ^o(r, 0)Ao(r', 0) + 

A colour-electric field strength is defined through 

g^Ei = i[DQ,Di]. (2.3) 
With this notation, the Euchdean correlation function considered [1] can be expressed as 

1 3-2. /ReTr \u{p,T)g^Ei{T,0)U{T,Q)g^Ei{Q,0) 



ReTY[[/(/3,0)] 



(2.4) 



where /? = 1 /T is the inverse temperature. Note that we keep a fixed "3" in the denominator 
even in the presence of dimensional regularization; this has no effect on the final renormalized 
result but does affect the finite parts of divergent intermediate expressions. 
The spectral function corresponding to G^{t) can be determined from 



Ge((x>„ 



/3 



dre*""^GE(r) , 



ImG'E(wn — ^ —i[uj + iO"*"]) , 



(2.5) 
(2.6) 



where the analytic continuation becomes unique by requiring powerlike (not exponential) 
growth at asymptotic frequencies; the momentum diffusion coefficient then follows from 



_ 2rpE(^) 

K = lim 

aJ->-0 OJ 



(2.7) 



If only a numerical determination of G^{t) is available, the task would be to invert the 
relation 



Ge(t) 







-da; 
— Pe(w) 
vr 



cosh I ^ — r I w 
sinh^ 



(2. 



for which various practical recipes have been proposed; see, e.g., refs. [TBj. (Probably the 
asymptotic large-w behaviour of Pe(w) is needed as input for the inversion.) 

The graphs contributing to the numerator of are shown up to 0{g'^) in fig. [TJ As 
is typically the case in finite-temperature QCD, however, certain subclasses of higher-order 
graphs may need to be considered as well. This can be achieved with the so-called Hard 
Thermal Loop (HTL) effective theory p!7]; where necessary, the observables will then be 
expressed as 
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Figure 1: The graphs contributing to the colour-electric field correlator, G'e(t) defined in eq. (I2.4p . 
up to 0{g'^). The big circle denotes a Wilson line wrapping around the time direction; the small 
dots the colour-electric field strengths; and the grey blob the 1-loop gauge field self-energy. Graphs 
obtained with trivial "reflections" from those shown have been omitted from the figure. 

where "naive" refers to an unresummed computation and "resum" to a resummed one. The 
unresummed difference in the first square brackets is infrared finite (provided that the correct 
low-energy effective description is used), and can be computed in naive perturbation theory, 
i.e. from the graphs of fig. [U in an effective field theory language, it corresponds to an 
ultraviolet matching coefficient. The last term of eq. ()2.9p is infrared sensitive, but can be 
addressed within the simplified effective description rather than with full QCD. 



3. Outline of the computation 
3.1. Leading order analysis 

In order to illustrate the steps of the computation, we start by explicitly deriving the (well- 
known) leading order result for the spectral function, corresponding to the contribution of 
graph (a) of fig. [U A direct evaluation of the Wick contractions leads to 

4'^(t) = ^ J]{a2A,,(r, 0) + dfAooir, 0) - dod, [Ao,{t, 0) + A,o{t, 0)] } , (3.1) 

i 

where A^^ is the gauge field propagator in configuration space. Expressing the propagator 
in momentum space, 

\uir,^) = i e*'="-+*'^-A^,(K) , (3.2) 
''k 

where ^j^- = TY^^,^ and is a bosonic Matsubara frequency, /c„ = 2TrTn, re G Z; and 
inserting the propagator of the general covariant gauge, 
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the longitudinal part proportional to the gauge parameter, ^, is seen to drop out. The 
remaining expression reads 



3 T 



K 



(3.4) 



where D = 4 — 2e is the spacetime dimensionality. Writing k'^ = — k^, omitting the 
scaleless integral 1 due to properties of dimensional regularization, and making use of 



ryikfl 7 



n^jk) 

kl+k'^~ 2k 



,W~r)k _^ gTfc 



, 0<T</3, 



(3.5) 



where k = |k|, we get 



glCF{D-2) 
6 



k n^{h) 



The Fourier transform as in eq. ()2.5p yields 

glCF{D-2) 



6 



k 



1 1 

+ 



(3.6) 



k — iujn k + ia;„ 



(3.7) 



and taking a subsequent discontinuity according to eq. (j2.6p we obtain 

glCF{D - 2)tt 



6 



k 



6{k -oj) - 5{k + u;) 



The final integral over k can now be trivially carried out. For D = A the result reads 

(2), . g'^Cp 3 , . 

Pe (w) = — — uj-\ (3.9) 



Gvr 



where we also expanded the bare gauge coupling according to eq. (j2.ip . 



3.2. Wick contractions at next-to-leading order 

Proceeding to the level 0{g^), the structures in the graphs (b)-(k) of fig. [T] need to be 
considered. (There is also a disconnected contribution originating from the numerator of 
eq. (j2.4p . as will be specified presently.) Given that the observable is gauge invariant, the full 
result must be independent of the gauge parameter ^; nevertheless the individual diagrams 
contain non-trivial terms proportional to ^ and i^^. We have checked explicitly that both 
structures cancel in the sum, and display intermediate expressions only for the Feynman 
gauge, = 0- (As a further cross-check, the small-frequency limit will be compared against 
that of an independent computation in the Coulomb gauge, cf. sec. 14. 3p . 
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Considering first the structures (b)-(f) of fig. [H and combining them with the contribution 
of the denominator of eq. (j2.4p , viz. 



ReTr[?7(/3,0)] 



1-glCp / dr/ dT'Aoo(r',0) +0{gi), 



(3.10) 



we obtain 
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{D - l)dl + G(t, 0) / dr' / dr"G(r' - r", 0) 



K 



L 



-{D - 1)^2 + {D- 2)/fc2 



Q 



(3.11) 



Here we denoted G{xo,x) = [Aqo{xo,x)]^^q, i.e 



G{xo,x) 



K 



and moved to momentum space, making use of 

'dr' e*«"^' 



A-2 I l<-2 



1 



(3.12) 



(3.13) 



etc, where qn is a bosonic Matsubara frequency and 5^,^ = (5n,o is a (periodic) Kronecker 
delta function. It can be verified that the Matsubara zero-mode part of eq. (I3.13P does not 
contribute in dimensional regularization in eq. (|3.1ip : the prime in indicates that it is 
omitted. The term with {D — 1)K'^ in eq. ()3.1ip does not contribute either, because the 
scaleless integral 1 vanishes in dimensional regularization. 

The contributions of the graphs (g) and (h) of fig. [T] can be worked out in a similar way. 
A straightforward computation yields 



giCpN.iD - 1) 



a^G(T,0) / dT'G(T',0) 



2giCFN^{D-l)^ ikne 



K 



(3.14) 



gtCpN^iD - 1) 



G(t,0) 



giCpN^iD - 1) 



K,Q 



Q^{K-Q) 



2 • 



(3.15) 



The graphs (i) and particularly (j) of fig. [T]are somewhat more complicated than the ones 
discussed so far. For (i) a few steps lead to the momentum space expression 

D-l (L»-2)fc2 



K ■^Q 



(3.16) 
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while for (j) it is perhaps simplest to remain in configuration space: 
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(3.17) 



Here we carried out a partial integration and made use of (9g + V2)G(xo,x) = —5^'^\x), in 
order to identify a term that cancels against that in eq. p.l4p : and expressed the remainder 
through the integral 



Mr) 



dr' 



dr' 



diG{xo - t',x) doG{xo - r, x)9jG(xo, x) - diG{xo - T,x)doG{xo,x) 
+ (D- 2)G(xo - r , x) \doG{xQ - r, x)V^G{xo, x) - V^G{xo - r, x)aoG(xo, x) 



(3.18) 



The subscript "5" refers to the fact that this is the most complicated among the 5 independent 
structures that appear in our result (the other ones are defined in eqs. (I325|)-(I32H|))0 A 
corresponding momentum space expression, albeit with no care taken of the Matsubara zero 
mode contribution, is given in eq. ()A.30p of appendix A. 

Inspecting finally the self-energy contribution, graph (k) in fig. [U we recall that making 
use of the substitution Q ^ K — Q in order to simplify the numerator, the Feynman gauge 
self-energy has the form 



glNc 



Q 



i. 



QHK-Qr 



{Q} 



QHK-qy 



and the corresponding contribution to Ge{t) reads 



(3.19) 



(3.20) 



Following the usual convention, ^{qj in eq. (|3.19p means that the Matsubara frequency is 
fermionic, i.e. = (2n + l)7rT, with n G Z; also, in eq. (13.201) . K = {kn, ki). 

Combining eqs. (j3.19p and (I3.20p . the transverse projectors eliminate the parts proportional 
to K^Ky from 11^;^, while the parts proportional to Qf^Qu can be expressed as 

'^{kn6i^ - ki5Q^){kn5iy - ki5ou)QfMQu = {K - Qfknqn + Q^kn{kn - Qn) + K^qniq-n " ^n) • 

(3.21) 



^In fact I5 could be subdivided into further independent structures, but here we keep it as one entity. 



7 



The first two terms are odd in the summation variable qn or fc„ — Qn, respectively, and give 
no contribution in eq. (j3.20p . so that only the last term matters. In addition there is a 
contribution from the term proportional to 6^^^ in U^i, which we write as 



(3.22) 



In total, then. 



2P e 

K 



iknT 



X 



3 Tl 
2 



+ 



4(D - 2)A;2 4.{D - 2)qn{qn - kn) 



+ 



Q^{K-QY ' K^Q^{K-Qy K^Q^{K - QY 
2{D-2)(D-l {D-2)k'^^ 



-N{ 



{Q} 



{D-l] 



K2 (1^2)2 

{D - 2)A:2 



Q^{K-Qf ' K^Q^{K-QY ~ K^Q^{K-Qf 



+ 



2 fD-l {D- 2)k^ 



(3.23) 



We now sum together eqs. (I3.1ip . (I3.14p - (|3.17p . (I3.23p . It can be noted that, in dimensional 
regularization, 



D-l 



1 



(3.24) 



,^yki + k^Y 2 kki + k^- 

As a consequence the factorized structures on the 3rd and 5th rows of eq. (I3.23P are propor- 
tional to 4 - D = 2e, and since the sum-integrals ^gl/Q^ = rVl2 + 0{e) and ^^pjl/Q^ = 
—T'^/24: + 0{e) are finite in dimensional regularization, there is no contribution from these 
terms. Defining furthermore 



Mr) 
Mr) 
Mr) 
Ur) 



K,Q 



K,Q 



Q^K-Q)^ ' 



ikn T 



K^Q^K-Q)^ ' 

K,Q ^ ^ ^' 
^2 ^hfii 



K 



(3.25) 
(3.26) 
(3.27) 
(3.28) 



and correspondingly 2^{2}i ^{3} ^oi the cases that Q is fermionic, and recalling X5 from 
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eq. (j3.18p . we can write the complete (unresummed) contribution as 



QCD, naive 



^ AD -2) y-Z2{T) - 2X3(t) +X4(r)J + -Z^{t) 



+ 



{D - l)X|i}(r) -{D- 2)X{2}(r) + 4X^3} (r) 



(3.29) 



3.3. Matsubara sums, spatial integrals, spectral function 



The next steps are to carry out the Matsubara sums over kn , Qn (generalizing eq. p.Sp ) ; 
to Fourier transform with respect to r (eq. ()2.5p ): to take the discontinuity across the real 
axis (eq. (|2.6p ): and to carry out the remaining spatial integrals. We illustrate these steps 
for one of the structures appearing, I2 of eq. p.26p . in some detail in appendices IA.lf[A.3[ 
An alternative derivation of the "vacuum part" of the most complicated structure, X5, can 
be found in appendix IA.4I The final results are collected in appendix IA.5I Inserting the 
expressions from eqs. ()A.54p - ()A.58P of appendix IA.51 together with D = 4 — 2e, we thus get 
the unresummed bare contribution to the colour-electric spectral function: 



QCD, naive 



+ 



2(1-6) 



-12{0J)-2I3{UJ)+I4UJ) +^X5(w) 



(3 - 2e)X{i}(u;) - 2(1 - e)X{2}(a;) + 4X^3} (, 



OJ 



(3.30) 



+ 



3 
1 

47^3 



-4€ 



(47r)3 



22/1 , 364 167r^ 

— - + 21n-p^ +— — 

3 V e 4a;2 / 9 3 



{q^ + 2uj^)ln 
1 



q + w 



w 



+ — F 



-it 



+ 



gtCpN^iu^ 
3 1 (47r) 
1 f°° 



4/1 

-3(7 + ^^^ 



In'- 


+ 0; 


+ 




u 




\ 52 


4a;2 


/ ~~9 



+ quj ^In ■ 
1 



UJ 



q — UJ 



In- 



q + w 


+ quj(^ 




q — w 





;( In- 



7^ — a;2| 



(3.31) 



To renormalize this, we need the generalization of eq. (13. 9p with effects of 0{e) included, 
which can be obtained from eq. (IA.2ip . taking into account the D — 2 from eq. (j3.8p : 



glc 



QCD, 



(3.32) 



Inserting here the bare gauge coupling from eq. ()2.ip and re-expanding in terms of the renor- 
malized gauge coupling, we get an additional contribution which removes the divergences 
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from eq. ()3.3ip : 



22Nc 



e 4a;2 



7^ + 1 • (3-33) 



QCD, naive 3 ' 3 

The sum of eqs. (j3.3ip and (j3.33p is our full result for the umesummed next-to-leading order 
contribution to Pe{<^)'-, it is reproduced in complete form in eq. (|4.2p below. 



3.4. Hard Thermal Loop resummation 

The result obtained so far (given explicitly in eq. (j4.2p below) contains no infrared divergences, 
so that in the regime uj>T there is no need for a resummation to be carried out. Nevertheless, 
we would like to extend the applicability of the result also to the regime uj < gT in which 
collective phenomena like Debye screening play an important role. The reason for this desire 
is on one hand formal, allowing us to cross-check that the framework is in principle well- 
defined in this limit as well; and on the other hand practical, permitting us to compare with 
existing results in the literature. In the regime uj < gT, logarithmic divergences do appear 
and a resummation becomes necessary. 

The way to consistently carry out the resummation is dictated by eq. (|2.9p . The practical 
computations are somewhat cumbersome; given that they also have no bearing on our main 
result (i.e. the behaviour of the spectral function in the regime uj>T), we relegate the 
discussion to appendix B. The main results can be found in eqs. (j4.3p . (|4.6p below, with a 
notation as introduced in eq. (|4.ip . 



4. Results 

In this section we collect together the results from the previous section and the appendices, 
and write down our final formulae for the colour-electric spectral function. 



4.1. Full expression 

In order to write down the full result, we rephrase eq. (j2.9p as 
Pb(w) = (Pe) -(Pe) + [P 

\ / HTL, resum 

+ (pe) + 

\ / IR, resum 



(pe 


) QCD 






(pe 


) QCD 






Pe) 






QCD, naive 



Pe 
Pe 

+ (Pe 



Pe ) - Pe 

HTL \ / IR. 



Pe ) + Pe 

HTL, naive \ / HTL, resum 



Pe 



IR., resum 



eq. (02 



eq. dO 



eq. dM 



(4.1) 



where we left out the arguments io for brevity. The term (pe)jp^ resum' ^^^ed and subtracted 
in eq. ()4.ip . is defined to be the HTL-resummed spectral function but computed only up to 
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the leading (linear) order in a Taylor expansion around u = 0. The rationale for this split-up 
is that all three structures are ultraviolet and infrared finite (i.e. contain no 1/e divergences) 
and that the last term, eq. (|4.6p . is readily available from the literature (cf. below). 
The first term of eq. (gT]) is the sum of eqs. ([MI]), ([3:3311 : 



Pe(w) 



QCD, naive 



1 



-N{ 



2 

— in — ^ 

3 iu^ 



2 20 
+ ^ 



2 149 
+ — 



3 



9'Cf 





27r2\ 






+ iVf^ 


/•oo 

/ dqnp{q) 





q + w 
q — w 




POO 

/ dqriBiq) 
'o 


(g2 + 2a;2) In 


q + w 
q — w 



+ qu!\ In 



+ — F 

q \q + UJ 



In- 



qoj yhv ■ 
1 



+ 



■In- 



UJ 



UJ 



■ UJ 



UJ\ 



(4.2) 



We stress that this spectral function contains no "transport peak" ~ uj5{uj) (in fact none 
of the master sum- integrals appearing in our computation have such terms). In addition, 
the function G^{ujn)-, the cut of which determines /Oe(ci;), contains no term constant in w„, 
which would correspond to a contact term in the Euclidean correlator Ge(t). Therefore p^{uj) 
contains all the information characterizing the correlator in an "unproblematic" form (see 
ref. [H] for a recent discussion of this issue in a context where complications do appear). 
The second term of eq. (|4.ip is a combination of eqs. ()B.6p . ()B.14p . ()B.20P : 



Pe(w) 



IR. rcsum 



HTL, naive 



67r 



where is the Debye mass parameter, 



— H 

6 3 



1 



wrrig I In 



2uj 

THe 



1 



(4.3) 



(4.4) 



The term of 0{g^) cancels in the sum of eqs. ()4.2p . ()4.3p : this term is contained in a modified 
form within eq. ()4.6p (cf. eq. (14. 7p ). It is also interesting to note that, inserting m| as 
expressed in eq. (IB.Sp . the next-to-leading order contribution in eq. (14. 3p can be rewritten as 

g^Cp g^ 

= — —r I iiq i\fUp{q)-tiycii'B{q) ywim 

HTL, naivo DTT ZTT^ ./fl L J V mt 



IR, rcsum 



which nicely combines with some of the terms in eq. ()4.2p 



/ 4^2 
quj m — — — 2 



(4.5) 
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The remaining ingredient, the HTL-resummed contribution beyond the hnear term in w, 
i.e. (/Oe)htl, rosum " (Pe)ir, rcsum, was determined in ref. [3]. With a minor change of notation, 
eq. (3.8) of ref. [4J can be expressed as 



HTL, rcsum 



IR, rcsum 



'dkk 



A;2 



^2 + 1 



(li^ I a) I 1 



, lni±^^ 



2 / . \ 2 



+ 



fc2 



+ 



+ 



+ 



1 



2fc fc-Lj / V 2fc 



(i|3(4 -(^2) + i| 



. In ±tiE = 



(4.6) 



'Ih p oj — fc I 



where w = uj/m-^ and /c = k/m^, and the expression within the curly brackets is dimensionless. 

For a) ^ 1 it can be checked, perhaps most easily numerically, that the expression within 
the curly brackets behaves as 



{...} 



(4.7) 



This cancels against eq. (j4.3p . so that for u ^ ttt-e only the naive QCD contribution of 
eq. (|4.2p is left over, as must be the case. 



4.2. Ultraviolet asymptotics 

An important application of the general result of sec. 14. H is that it determines the asymptotic 
(oj ^ {T, mE,At^}) behaviour of the colour-electric spectral function. In fact, as explained 
in connection with eq. (|4.7p . for large frequencies to S> 171^ we can make directly use of the 
unresummed expression, eq. (j4.2p . Expanding the integrand of eq. ()4.2p in q/co and making 
use of 

120 ' 



6 



1 



vr 
15 



+ 1 



(4.^ 



it can be seen that the thermal part of the spectral function disappears at large frequencies: 



Pqcd(w) 



T=0 



+ 



67r 180a; 



(7iVf - lliVe) . 



(4.9) 
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Such a power-law decay of the thermal correction at ^ T is in accordance with the general 
results of ref. [19], although it is not clear to us whether a precise relation to the Operator 
Product Expansion as used in ref. [19] can be established for our non-local correlator. The 
term denoted by "T = 0" in eq. ()4.9p . on the other hand, is temperature independent, and 
given by the first two rows of eq. (|4.2p : 



T=0 



1 + 



9 



(4vr) 



2 p,^ 20 



, 11 p^ 149 

+ A^c — In H 



3 



+ . (4.10) 



A rough estimate can be obtained by setting fj, ^ 2uj and inserting the 1-loop running coupling. 



Pqcd(w) 



47rCf 



T=o (llA^c - 2iVf) ln(2a;/A^ 



1 



(247r2 - 149) TVc + 20iVf 
6(lliVc -27Vf)ln(2tj/AMs) 



(4.11) 



Numerical evaluations are discussed in sec. 14. 4[ 

It is perhaps appropriate to stress that the first omitted term in eq. (|4.7p (which is not ex- 
plicitly cancelled within our computation) leads to a contribution of the type dp^ ~ g'^m'^/uj, 
which is parametrically of higher order than the T-dependent part of eq. (j4.9p . Hence, at 
0{ag), ultraviolet asymptotics is indeed completely determined by hard modes. This is again 
in accordance with the general analysis of ref. |19j 



4.3. Infrared asymptotics 



We next discuss the infrared {uj <^ ttie) behaviour of the spectral function, in the range 
T ^ Aj^. Although we have nothing to add to the existing results in this regime [TJ [8] , it 
is comforting to crosscheck that we can at least reproduce the expression at 0(0^). 

bmce 

the computation of ref. [7| was carried out in Coulomb gauge whereas we have described a 
covariant gauge computation, this also serves as a further cross-check of gauge independence. 

In the regime w <C ttie the "T = 0" contribution (eq. (|4.10p ) is subdominant compared 
with the thermal one, which is of O^g^m^uj). Expanding the thermal integrand of eq. ()4.2p 
in uj/q and making use of 



dx X 


7r2 


dx X In X 


_ 7r2 


- 1 


~ 6 ' 


io - 1 


6 


dx X 


_ 7r2 


dx X In X 


_ 7r2 


+ 1 


~ 12 ' 


Jo + 1 


~ 12 



1-7E + 



m 

C(2) 



In 2 + 1 - 7e + 



m 

C(2) 



(4.12) 
(4.13) 



we obtain the infrared behaviour of the unresummed QCD result, expressed now in terms of 
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the momentum diffusion coefficient of eq. (|2.7p : 



]qCD, naive — 



u 

T 3 

ln- + - 
m 2 



QCD, naive 

C'(2) 



7e + 



C(2) 



■ 2T 3 

In— + 
w 2 



7e + 



C(2) 



• (4.14) 



Tlie symbol "~" signals that the logarithmically divergent term has been kept finite on the 
right-hand side. Similarly, from eq. (j4.3p . the HTL contribution to k is 



JHTL, resum 



J HTL. naive 



(4.15) 



Summing together eqs. (|4.14|) and (j4.15|) . we reproduce the result of ref. [7]: 
g^CpT of 2T 1 C'(2) iVfln2 \/ , 



(4.16) 



As indicated, corrections to this expression start already at 0{g), and have been determined 
in ref. [8], but in the infrared regime these go beyond the accuracy of our computation. 

We note that if g is not small {T/m^ is not large), the expression in eq. (I4.16P can become 
negative. This unphysical behaviour is a result of extrapolating the asymptotic weak-coupling 
expression beyond its range of applicability; for A'c = 3, A^f = 0, the problem occurs already at 
g > 1.05 and is another reflection of the exceptionally poor convergence of the weak-coupling 
expansion of k (see also ref. |8j). Given that the unresummed expression (eq. ()4.14p ) does 
stay positive and in fact diverges as a; — )• 0, the HTL-resummation (eq. (|4.15p ) in some sense 
"over-corrects" the result in this regime. 



4.4. Numerical evaluation 

In order to finally evaluate our result (eq. (14. 2j) -|- (14. 3j) -|- (14. 6p . cf. eq. (14. ip ) numerically, we 
need to insert some value for the running coupling g^. It is only in the regime uj ^ T that 
two subsequent orders with the same functional form are at our disposal, so that some kind 
of a scale optimization is possible; we can then define /iopt for g^ from the "fastest apparent 
convergence" or "principal of minimal sensitivity" criterion based on eq. (|4.10p : 

lnfc,,„,).ln(2.) + ^^^^^j^^^i-^^^, (4.17) 

a structure indirectly also visible in eq. (j4.1ip . In the infrared regime oj <^T such an exercise 
is not possible but we nevertheless need a value for g^; we choose it from a context in which 
a next-to-leading order computation does exist, namely that of "EQCD" (cf. ref. pO] and 
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Figure 2: A numerical evaluation of p-e{^), eq. (|4.ip . in units of wT^, for Nc = 3,Nf = 0, T = 
S.TSAms (left) and T = ISA^g- (right), corresponding to T « 3Tc and T « 12Tc, respectively. The 
gauge coupling has been fixed as explained around eqs. (|4.17l) . (|4.18p . and the error bar reflects the 
corresponding uncertainty. (In the 0{g^) result the "optimal" scale is fixed to the thermal value of 
eq. (|4.18p . i.e. does not change with the frequency.) A slightly negative intercept at w — > is an 
artifact of the truncated weak-coupling expansion, cf. discussion below eq. (|4.16p . 



references 



therein) |§ 



ln(^opt(T)) = ln(47rT) - 7e 



Nc-8 In 2 iVf 



(4.18) 



2(lliVc - 2iVf) 

Comparing eqs. (j4.17p and (j4.18p we see that we can switch from one regime to another at 
a; ~ T, for instance at a; ~ 0.8903T for Nf = 0. In order to get an error band for the 
uncertainty related to the choice of g'^, we use the 2-loop g'^, with the scale parameter varied 
in the range (0.5 . . . 2.0) x fiopt- Since the Debye mass parameter plays no role in the range 
CO ^ T that we are mostly interested in, we fix it in the simplest way imaginable, through 
eq. (|4.4p with the gauge coupling evaluated as described above. 

The outcome of a numerical evaluation of eq. ()4.ip is plotted in fig. O in units of wT^ . We 
observe that because of the running of the gauge coupling the next-to-leading order result 
falls below the leading-order one in the ultraviolet domain, u ^ T. The dependence on the 
scale choice is also drastically reduced, becoming (perhaps rather surprisingly) practically 
invisible as soon as uj>T. 

In the infrared domain, <C T, the next-to-leading order correction eventually overtakes 
the leading-order result, and even the HTL-resummed expression shows very slow convergence 



^For EQCD even a NNLO computation exists [21] but that goes beyond the accuracy of the present study. 
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Figure 3: A numerical evaluation of Ge(t), in units of T'^, for = 3,Nf = 0, T = 3.75A^js (left) 
and T = 15A^ (right), corresponding to T « 3Tc and T w 12Tc, respectively. The gauge coupling 
has been fixed as explained around eqs. (j4.17p . (I4.18|) . and the error bar reflects the corresponding 
uncertainty. (In the 0{g^) result the "optimal" scale is fixed to the thermal value of eq. (I4.18p . i.e. 
does not change with the frequency.) On the logarithmic scale used a clear temperature-dependence 
is visible only around the middle of the interval, r T « 0.5. 



(cf. ref. [8j). Nevertheless, fig. [2] does illustrate an important point: even at weak coupling 
there is no transport peak around the origin; rather p^{uj)/uj displays a relatively flat behavior 
at CO <T. More precise studies of the infrared domain have so far been carried out within 
classical lattice gauge theory [11] as well as within = 4 Super- Yang-Mills theory at infinite 
't Hooft coupling [3 [TO]. 

In fig. [3 we plot the Euclidean correlators G'e(t) corresponding to the spectral function 
Pe{^jj) of fig. [H obtained through eq. ()2.8p . The leading-order result can be given in a closed 
form |4], 



g CpTT I 



+ 



1 



cos^ (vttT) 

sin^(7rrT) ' 3sin'^(7rTr) 

(4) 



(4.19) 



whereas the next-to-leading order correction has been evaluated numerically. The ui- 
dependent difference of at 0{g'^) and 0{g^) in fig. [2] has converted to a practically r- 
independent shift on the logarithmic scale of fig. [3| 

We note that, parametrically, the range uj < <^ 2T gives a contribution to eq. (j2.8p of 
magnitude 

f""^ da; TpeH 
X 70 vr oj 



g^mlT ~ 0(<7°r 



(4.20) 
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cf. eq. (j4.15p . which is formany of higher order than the ultraviolet contribution, 0{g^T'^). 
This is also reflected through the fact the infrared divergence of eq. (|4.14p is integrable at small 
uj. Nevertheless, the next ultraviolet contributions are of higher order still, 0{g^T^), so we can 
already meaningfully include even the infrared part of Psi^) in our numerical evaluation of G^- 
In practice, of course, the scale ttie is not smaller than ~ 2T at phenomenologically interesting 
temperatures, and the infrared contribution should not be substantially suppressed. 



5. Colour-electric correlator in the strong-coupling expansion 

In the confined (low-temperature) phase of pure SU(A''c) gauge theory, the Euclidean corre- 
lator of eq. (|2.4p looks singular, because the expectation value in the denominator vanishes. 
We wish to demonstrate that despite this appearance the ratio could actually remain finite, 
if defined through a suitable limiting procedure by explicitly breaking the Z(Nc) center sym- 
metry (an obvious possibility would be the addition of dynamical quarks). As a tool we use 
the lattice strong-coupling expansion, which converges (deep enough) in the confined phase. 
The result also yields an interesting functional dependence on r, which can be contrasted 
with the weak-coupling result of fig. [3l 

Defining the strong-coupling expansion requires the use of lattice regularization, and hence 
also a discretization of the correlator of eq. (j2.4p . The discretization is obviously not unique; 
we choose the possibility proposed in eq. (4.2) of ref. [1], viz. 

E L 1 Re ( _ ( r-j ) ~ ( L-n ) +~ ( L-n ) _ ( r- J ) > 

G'E(r) = xo 

-6a^ReTr(. ■) 

(5.1) 

where o is the lattice spacing and we have used a graphical notation with lines indicating 
parallel transporters (for more details see ref. [3]). Let us also define the auxiliary correlator 

C(t) ^ReTr( | \) 

\ (5.2) 

/3-r r 

in which the vertical lines have the length of one lattice spacing. If C(0) 7^ 0, we can then 

write . . , . 

r I ^ C{T-a) + C{r + a)-2C{T) 

G.ir) = . (5.3) 

Note that this expression makes sense only for t > a. 

In order to guarantee that C(0) is non-zero, we add a small "source term" to the (Wilson) 
action, thus breaking the Z(A'c) symmetry explicitly. We believe that the precise from of 
this regulator does not matter, and in practice we choose the simplest possibility, adding to 
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the action just the volume sum of the Polyakov loop itself, multiplied by a small "magnetic 
field", to be denoted by k. Then (if iV^ > 4), 

C{t) = K u-a {l + 0{u^)^ +0{k'^) , (5.4) 

where A',- = /3/a and u is a variable naturally arising in the character expansion (for details 
see, e.g., ref. p2]). 

uEEi-Tr((P,,)). (5.5) 

Here Pij denotes an elementary plaquette in the (i,j)-plane, and the expectation value is 
defined as {{Pij)) = J dP^, exp(|^TV [P^, +4])/ / dP^,- exp(|^Tr [Pi,+4.]). The strong- 
coupling expansion is an expansion in a small Wilson parameter /3g, which means that the 
variable u should also be small: for A'c > 3, 

+0{/3l). (5.6) 



2iV2 



Incidentally, we have also worked out the correction of 0{u'^) to eq. ()5.4p : it amounts to the 
replacement — )• (1 + 4:xu'^)u^, with x = ^ and x = but this does not change the 

qualitative behaviour so we do not elaborate on further details here. 

Inserting eq. (j5.4p into eq. (|5.3p . we note that the magnetic field k drops out, and a finite 
result is left over. Moreover, inserting some small value for w (0 < n <C 1), we obtain 
a function which is not unlike those in fig. [Sj examples are shown in fig. [H The main 
differences are that at very low temperatures (very small u), the exponential suppression 
around the middle of the interval is much stronger than in the deconfined phase; and that 
there is little curvature in the plots (i.e. only a single exponential). Taking a discrete Fourier 
transform and subsequently naively applying eq. (12. 6p to it, leads to a spectral function with 
(5- function peaks at w = it^ Inu but no other structure. This behaviour is completely different 
from those in fig. [21 the difference demonstrates once again that a high resolution is needed 
in the determination of Gb(t), in order for us to be able to extract the correct non-trivial 
features of the corresponding spectral function. 

In a practical lattice measurement, the numerator and the denominator of eq. (j5.3p are 
evaluated separately. If the breaking of the Z{Nc) symmetry is only explicit (not sponta- 
neous) and small, the signal is likely to be very noisy. This is not really a problem from 
the phenomenological point of view, though, since heavy quark diffusion appears anyways 
to be a physically meaningful concept only in the deconfined phase where 7j{Nc) is indeed 
spontaneously broken (at any non-zero lattice spacing). 
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Figure 4: A numerical evaluation of the strong-coupling result for Ge{t), eqs. (|5.3p . (|5.4p . in units of 
T^. For a fixed A^r, smaller values of u correspond to lower temperatures. The qualitative behaviour 
can be compared with that in fig. [3l the exponential suppression around the middle of the interval 
becomes very strong at low temperatures (note that the scales are different in the two figures). 



6. Conclusions and outlook 

The "transport coefficient" related to the correlator of two colour-electric fields along a 
Polyakov loop determines the momentum diffusion coefficient, n, of a heavy quark near rest 
with respect to a thermal medium. Among all the QCD transport coefficients this is probably 
the most accessible one from the theoretical point of view; indeed it is the first for which a 
next-to- leading order weak-coupling expression was obtained [S]. Alas, as the weak-coupling 
expression itself ^ as well as numerical studies within classical lattice gauge theory [11] show, 
the weak-coupling expansion converges very slowly at any temperature relevant for heavy ion 
collision experiments, so a non-perturbative determination would be highly desirable. 

On the point of the non-perturbative determination, there is the salient feature that no 
transport peak exists in the colour-electric spectral function from which k is extracted, /Oe(w), 
even at weak coupling where spectral functions are in general more singular than in a strongly 
interacting system. The reason is that there is no (approximately) conserved charge related 
to the colour-electric field strength. Employing this theoretical insight as a constraint might 
help to reduce the systematic uncertainties that are inherent to the practical recipes used for 
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the analytic continuation from the measured Euchdean correlator to the desired Psi^)- 

On the other hand, the non-perturbative study is faced with the challenge that the absolute 
magnitude of the correlator is needed for the determination of k, which is obtained from the 
a; —7- intercept of 2Tp^{Lo)/bj. This means that the renormalization of the colour-electric 
field entering the correlator plays an important role. 

In principle, the colour-electric field can be renormalized independently of the particu- 
lar observable considered (perhaps in analogy with procedure for colour-magnetic fields in 
ref. [23]). However, its discretization is not unique, so it may be useful to have a way to 
"crosscheck" the normalization of the result directly with the data at hand. For this pur- 
pose our asymptotic large-w behaviour might turn out to be useful, since the weak-coupling 
result is relatively accurate at a; ^ T (at least within ~ 10 %, which would at the current 
stage be perfectly sufficient as far as k is concerned). Our formula is strictly applicable only 
in the continuum limit; nevertheless, a practical matching could already be attempted at 
some intermediate frequency a; and non-zero lattice spacing a, provided that the inequali- 
ties T ^ w <C 1/a can be satisfied. Of course, all of this assumes that some p^{ijo) can be 
approximately reconstructed from a measured Euclidean correlator Ge(t) [16j. 

Apart from the spectral function (fig. [2|), we have also estimated the corresponding Eu- 
clidean correlator (figs. El [1| which can be directly compared with lattice simulations. The 
result appears to be most sensitive to thermal effects around the middle of the Euclidean time 
interval, where its absolute value is very small, posing a challenge for numerical simulations. 
Nevertheless, with refined techniques a reasonable signal should ultimately be obtainable. 

To summarize, it is our (perhaps optimistic) hope that by combining our analytic work with 
forthcoming numerical studies a qualitative and even approximate quantitative understanding 
concerning the shape of the spectral function p^{ijj) and in particular the intercept n = 
lim^^^^o '2Tp^{uj)/uj can eventually be developed. 
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Appendix A. Details of the unresummed computation 



A.l. Matsubara sums, Fourier transform, spectral function 

We illustrate here the steps needed for the determination of the spectral function with the help 
of one of the basic structures appearing in the analysis, the function defined in eq. (j3.26p . 
It turns out that in the course of the computation the expression splits into several parts, 
some of which may contain poles for specific spatial momenta; although such poles cancel at 
the end, it is important to "regulate" the intermediate expressions consistently so that no 
finite pieces go amiss. A convenient way to accomplish this is to redefine as 

K,Q 

We have considered another regularization as well, writing the denominator as 1/{K'^[Q'^ + 
X'^][{K — Q)"^ + A^]}, and checked that this leads to the same result in the limit A — )■ 0. 
As a first step, we carry out the Matsubara sums over A;„,g„. Defining 



Ek = + A2 , Eq=H, Efcg=|k-q|, (A.2) 
the sums can be written as 

{kl+El){ql + El)[{K-qr:)^+El^] 

^kl+El){ql + El){rl + El^) 

/3 ^ gifcn(T+o-) ^—iqncr ^ir„a 

where in the last step we used a representation of the (periodic) Kronecker delta, /36t^ = 
J^dae'^^". The sums have now factorized and can be carried out like in eq. p.Sp : care must 
be taken just in order to transform the time arguments to the interval (0,/?), which can be 
achieved for 

qn 

by making use of the symmetries 

G{-T,E) = G{T + mf5,E) = G{T,E) , mGZ. (A.5) 

Since r, a only appear inside exponential functions in the result (cf. eq. (|3.5p ). the integral 
over a as well as the subsequent Fourier transform (eq. ()2.5p ) are trivial. The discontinuity 
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as in eq. (j2.6p can finally be taken by first setting exp{ziziu}nP) = 1 in order to remove terms 
growing exponentially at large |a;|, and by subsequently making use of 

1 



Im 



(A.6) 



Note that by making use of this relation we are advised to replace all other fractions by 
principal values {t^^^ = F(i) - i7r(5(x)). 

Let us stress that all terms (in this as well as in other master sum-integrals) contain 
exclusively within the structures of eq. (lA.Gh . always with at least one non-zero Ei; therefore 
no "transport peaks" ~ u;6{uj) arise in the analytic continuation. In addition, there are no 
a;„-independent terms, which would correspond to contact terms in configuration space. 

Denoting the result of these steps by X2(a;), we obtain 

TTk^ 



T2OJ) 



lim 



Jk.q ^EqEf^q 



F 



^ 6{uj - Ek) - 6{uj + Ek) ^ 



5{uj - Eg- Ekq){l + riBi + nea) 

+ 5{uj + Eg - Ekq){nBi - nB2) 
+ 6{uj - Eg + Ekg){ns2 - nsi) 

- 5{uj + Eg + Ekq){l + nei + n^2] 

1 1 



2Ek 



2Eg n^2 



Ek + Eg + Ek 
1 



+ 



-Ek + Eq + Ek 
1 



{Ek 



E,? 



Elq 



+ 



{Ek+EqY-E^ 



+ 



{Ek-EkqY-El {Ek + EkqY-El 



where 



n^{E^ 



Up 



n^{Ekq) ■ 



(A.7) 



(A.8) 



The expression is obviously antisymmetric in w — )• — w; in the following we restrict, without 
loss of generality, to a; > 0. 

In eq. ()A.7p . two different qualitative structures can be observed. We refer to the first 
group, with the constraints 5{uj — Eg — Ekq) etc, as "phase space integrals"; and to the second 
group, with the constraint 6{uj — Ek), as "factorized integrals". We now proceed to discuss 
these two groups separately. The analysis parallels that in appendix A of ref. 



A. 2. Phase space integrals 

The phase space integrals are ultraviolet finite so that we can set Z? — )• 4 from the outset. 
We fix the z-axis in the direction of k and carry out the integral over the directions of q by 
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dqq 



changing integration variables: 

f d^q 

J (27r)3 ^ ilT^k 

For oj > 0, tlie following results are easily verified: 

rk+q 



k+q 



\k-q\ 
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k—q 
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dEkq6{uj -Eg- Ekq) 4>{Ekq) 



{uj-k)e 



k-q\ 
LO + k 



dEkq5{0J + Eq- Ekq)(t>{Ekq) = e{k-u)e[q 



i-K+q 

/ dEkq 6{0J -Eq + Ekq) ^{Ekq) 
J \k—q\ 



e{k-u)i 



k — uj 
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k + LO 



(f){q + uj) 
4>{q - uj) 



(A.9) 

(w - g) , 
(A.IO) 
(A.ll) 

(A.12) 



where we also simplified the notation as Eg -+ q. 

To go further let us start by considering the temperature-independent "vacuum part" , to 
be denoted with 5psfil2- 



'5ps,oi'2(w) 



lim 



d'^k 



rP 



A^o J (27r)3^ \El-ij^^ 



TTk' 



d3q 6{iO -Eq- Ekq) 



(27r)2 



^EqEkq 



(A.13) 



Changing integration variables as in eq. ()A.9p and making use of eq. (jA.lOp , we get 



(^ps,o2^2(w) 



47r 
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A^o (27r)3 

lim ^ 

A^o 327r3 

lim 



dfcP 



dA;F 



7r/c3 
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fc2 + A2 - a;2 y (47r)2 
k^ 



Aq 



fc2 + A2 - 



a; 



, A2 8 
^"4^ + 3 



TTo (47r)3 V^^^4l^2 ' 3y • (^-14) 

We see how the phase space integration contains a divergence, originating from an "on- 
shell" configuration (A; ~ w); this will be cancelled against a corresponding divergence in the 
factorized vacuum integral, eq. (IA.25j) . such that in the complete result the limit A — t- can 
be taken. 

Considering the thermal part, the substitution q — t- k — q, i.e. Eq -H- Ekq, allows us to 
remove some redundancy, and the integral becomes 



^p^,t'I2{<-^) = lim 



A^'o (27r)3 

k+q 



dk¥ 



TTk^ 
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4^ 



dq 
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/ dEkq \25{uj -q - Ekq)n^{q) + 25{u + q- Ekq)[ns{q) - nB{Ekq)]\ 
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A^o 16tt'^ Jq 
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(A.15) 
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where we made use of eqs. ()A.9p - ()A.lip . and subsequently shifted q + u ^ q in one of the 
integrals, in order to always have q as the argument of the Bose distribution. At this point it 
is advantageous to change the order of integration; the integration area can be illustrated as 



CO 



Thereby we obtain 

(5ps,Ti'2(^^) 



1 



lim „ , 

A-s>o IGtt-^ Jo 

1 

16^. /o 



D /•L0+2q 

dqu^iq) / dkP 
dqrijiiq) 



Aqoj + — In 



A;2 + A2 - w2 
q + oj 



OJ 



(A.16) 



A. 3. Factorized spatial integrals 



The factorized integrals are more delicate than the phase space ones in the sense that the 
q- integral can contain ultraviolet divergences. In general this means that the integration 
measure is to be written as 



d'^q 



dq q 



d-l 



+1 



dz(l - 



(A.17) 
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where d = D — 1 and z = k • q^/kq. If the integrand is independent of z, the z-integral yields 

^^^^^^ 



dz(l - z' 



r(i)r(^; 



r(f) 



If the q-integral is finite, as is always the case with thermal corrections, it is convenient to 
change integration variables like in eq. ()A.9p . 



+1 2 ff'+Q 

dz ^ — dEkq Ekq . 

-1 Kq 



(A.19) 



Sometimes one is also in the lucky position that the q-integral can be directly identified with 
a known vacuum integral (see below). 

Let us start by considering the temperature-independent "vacuum part" , to be denoted by 

d'^k TTk'^6{uj - Ek) 



lim , , 

A^o J (27r)« 

d^'q 
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2E. 
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{2t:Y AEgEk, 
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-UJ + Ea + E, 



kq 
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Here we again restricted to u; > 0. Within the k-integral the hmit A — )• can be immediately 
taken, so that Ek — )■ k; making use of eqs. ()A.17p . (jA.lSp . we get 



TT 



d'^k k5{uj-k) 



3,,-2e 



l + dln£,+2 
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As far as the q-integral is concerned, we recall that one of the basic 1-loop integrals can, after 
integration over q^, be written as 

d^Q 1 



i?o(i^';0,0) 



{2'k)D Q2{K - QY 
d'^q 1 



Therefore the 2nd row of eq. (1X20]) equals 
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iko + Eq + Ekq 
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-iko + Eq + Ekq 



ReSo(-w^ + A;2;0,0) = Re Bo(-A^; 0, 0) = Re 



2.1 



-2e 



i + ln4, + 2 



(A.22) 
(A.23) 

(A.24) 



where we made use of 5{u: — E^) in eq. ()A.20p in order to write uj"^ = k"^ + A^. In total, then, 
6,M^) = lim f - + In ^ + In ^ + 4V (A.25) 



A^o (47r)2 



4tj2 



A2 



Here we observe again the significance of the intermediate regulator. 

Summing together eqs. (|A.14p . ()A.25p . logarithms of A cancel against each other, and we 
obtain the final vacuum part of X2, given on the first row of eq. (|A.55p . 

As far as the thermal part is concerned, we change integration variables on the last row 
of eq. ()A.7p so that the argument of is always Eq. Furthermore, since the integral is 
exponentially convergent, we can use the measure in eq. ([Al9]) . Partial fract ioning the 
integrand and inserting E^ ^ oj, Eq ^ q, we get 



lim , , 

k+q 

dEkqV 

k — q\ 



°° T:k'^5{uj - Ek) 2tt 
ak - 
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(27r)3 
1 



dqrijiiq) 
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(A.26) 



The integral over Ekq is elementary, and sending A — t- 0, such that uj = + A^ — )• k, it 
yields 

[{u;-qf-{k + qn{u + qf-{k + qY] 



lim In 

A^o 



In total, then. 
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q — UJ 



327r3 Jq 
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q — UJ 



(A.27) 



(A.28) 



Summing together eqs. (|A.16p . (jA.28p . we obtain the thermal part (second row) of eq. (|A.55p . 
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A. 4. An alternative derivation of the vacuum part 



While the method of sections IA.mA.3l works in all cases, the temperature-independent "vac- 
uum parts" of the results can also be derived in a somewhat more straightforward manner. 
We illustrate this for the most complicated master integral appearing in our study, I^ir), 
defined in eq. (j3.18p . The topologies from which this structure originates, already sketched 
in graph (j) of fig. [H can be depicted somewhat more specifically as 

r T 

(A.29) 

r I 



The idea of the approach is the following. Going to momentum space; carrying out several 
changes of integration variables; and ignoring possible problems with Matsubara zero modes, 
Qn = (this will be justified presently), the integrals over x and r' can be carried out, and 
we end up with 





Mr) 



K 

2kn 



Q2(^_Q)2 

Q2 + (D_2)(fc2 - 



{2D - 5)k^ + (-31) + b)k^ 



(A.30) 



We now replace the sum- integral through its zero-temperature limit, the D-dimensional 
integral Jq. In addition, the Fourier transform over r and the subsequent analytic contin- 
uation and cut amount, due to the symmetry of the expression in /c„ — t- —k^, to simply 
setting kn — )• —i{uj + iO'^) and taking the imaginary part. Thereby any dependence on the 
temperature disappears. Before the analytic continuation, it is useful to write k^ = — k'^. 

The only complication is the handling of the second row of eq. (IA.30p . Amusingly, the 
structure is similar to that met in connection with the so-called cyclic Wilson loop in ref . [25] If| 
We note, first of all, that in the zero-temperature limit the Matsubara zero mode plays no 
special role, and the apparent divergence l/q-n can be regulated by taking a principal value. 
Furthermore, introducing a Feynman parameter, the integrand can be written as 



2kr, 



Qn {kn - Qnf + (k - q)2 



kn 



1 



{k„ 



)2 + (k-q)2 



(A:„ + 9n)2 + (k - q)2 



[{kn - qn? + (k - q)2][(A:„ + qnf + (k - q)2] 
{(k - q)2 + [{2s - l)kn - qn? + 4s(l - s)kiy 



(A.31) 



^In that case the analysis was originally carried out by M. Vepsalainen. 
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The discontinuity corresponding to eq. ()A.30P can now be expressed as 
M.) = 2Im|^{/-^5j^J-^[5(D-2)K= + (-3I, + 5)A-= 



-1 

u2 



+4K / ds 



1 



2 



Q(Q2 + Af2) 

-4(1) - 2)klikl - k2) [\s [ —— i —2 [ , (A.32) 



where = 4^(1 _ ^^{.2 ^2 = - ifkl + k^. 
The Q-integrals in eq. ()A.32p are all familiar: 



^-2^ /I 



1 /I , 



- + + 2 ' (A.33) 



+ ln^ , (A.34) 



In ^- . (A.35) 



QQ2[(^_g)2 + ^2]2 (4^)2i^2 ^2 

In the middle case, the remaining integration over the Feynman parameter is also straight- 
forward: 

It remains to take the discontinuity and to carry out the integral over k as well as, where 
present, over s. For the structures on the first row of eq. ()A.32p this leads directly to the 
"known" cases in eqs. ()A.54p . ()A.55p . We discuss explicitly only the 2nd and 3rd rows of 
eq. (|A.32p . since these structures do not appear elsewhere. 

Like in sections IA.21 IA.31 poles may appear in individual terms but are cancelled at the 
end; to handle these, we regulate the expression by replacing — )• + and denoting 
El = k'^ + \^. Inserting eqs. (|A.34p . (IA.36p . the structure on the 2nd row of eq. ()A.32p then 
becomes (omitting numerical coefficients) 



We note that, for a; > 0, 



a;2 „ / \ ivro; 



- + ln- ^ + 2 = - + ln-p^ + 2 + i^. A.39 

e 4(a; + zO+)(-a; - ^0+) e 4l^2 ^ ' 
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A non-zero contribution only arises by combining the imaginary part of eq. (|A.38P with the 
real part of eq. ()A.39p . because the other possibility leads to 



F 



1 



P + A2 - a;2 



d'^k 1 



P 



A2 



47r2 



dkF 



1 



_ A2 k + Voj^- A2 



. 



(A.40) 



Making use of eq. ()A.2ip in order to evaluate the k-integral then leads to 

-,2 



(^2X5 (tj) 



(47r) 



(A.41) 



It remains to consider the structure on the 3rd row of eq. (]A.32p . Inserting eq. ()A.35P and 
shifting the Feynman parameter by | in order to make the integrand symmetric, it can be 
expressed as 



53X5 (w) 
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-2e 



(4vr) 



ds / Im 



k„ 



A-2 _ F 



fe2 + El {2knsY + A;2 



In 



ki 



ln(l - 4s^ 



. (A.42) 



The integrand is the product of three terms; the first one is like in eq. (IA.38p . while the 2nd 
and 3rd terms can for > be written as 



UJ^ + k^ 



-2ojs + k - i^+){2ijjs + k) 



F 



a;2 + A;2 
/c2 — 4a;2s2 



+ i'Ke{uj - k)6{2ujs - k) 



2k 



(A.43) 



ln(-" + ^T,-°")^" + ^^-ln(l-4.2) 
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■ UJ 



ln(l - 4s2) + inO^k - w) . (A.44) 



{uj + iO+){-uj - iO+) UJ- 
Various possible contributions can now be identified, according to whether the imaginary 
part comes from eq. ()A.38p . ()A.43p . or ()A.44p (there is no contribution from multiplying the 
three imaginary parts, because of the appearances of 9{uj — k) and 9{k — u) in eqs. (lA.43p . 
()A.44p . respectively). Actually, the imaginary part of eq. (fX43]) gives no contribution, be- 
cause the real part of eq. (IA.44P vanishes for the value of s that is allowed by the (5-function 
of eq. (IA.43p . This then leaves two possibilities. 

If we take the imaginary part from eq. (IA.44P and write 9{k — oj) = 1 — 9{uj — k), then the 
integral over unity can be reduced to integrals like that in eq. (fXiOj) and gives no contribution. 
A non-zero contribution comes from the integral (since the integral is UV-finite we set d — )• 3) 
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1 + X 



(A.45) 
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where we carried out the integration over s and denoted x = k/uj; the two terms inside the 
square brackets come from the partial fractions of the first principal value (we have set A — >■ 
wherever the limit exists). 

If we take the imaginary part from eq. ()A.38p . the integral becomes 



■KUJ 



(47r) 



ds 





'dk k"' 
'2^ 



|A;2_,.,2| 



K-^-Ek)^[ J T.:.„2 ) In^— 5^-ln(l-4s2) 

(A.46) 

This time the integral over k is trivial. Denoting x = 2s, a part of the x-integral can also be 
carried out, and we arrive at 
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647r3 \ 4tj2 
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+ F 
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v/l-AVw^-x 



ln(l - x^ 



(A.47) 

Adding up eqs. ()A.45p . ()A.47p . the coefficient of ln(l + x) develops a prefactor x"^ — 1; this 
cancels the pole (for \/oj — t- 0) and leaves a regular finite integral. The coefficient of ln(l — x), 
in contrast, has a prefactor + 3 = 4 + 2(x — 1) + (x — 1)^, so a singular integral remains. 
Denoting y = \/\ — y? jup- — x and 5=1— \/l — A^/a;^ ~ jluP' , the singular integral can 
be written as 
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idy 
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in h in y 



(A.48) 



where in the first integral of the last row we wrote x = y/6 and in the second could set 
(5 — )• in the upper limit of the integration. Of the remaining integrals the singular one is 
now trivially integrated. 



4/'|l„, = 21„=(*,.21„^(^) 



(A.49) 



whereas it requires a few more steps to show that 
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Summing up, all terms containing logarithms of A nicely cancel against each other, and we 
obtain 
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We note that the pecuhar vr^-terms, which appear nowhere else in our computation, can be 
traced back to the first term on the last row of eq. (|A.48P as well as the term in eq. ()A.50P : 



Q X 1 — X 3 



(A.52) 



Inserting finally the vacuum parts of eqs. (|A.54p . ()A.55P as well as the newly determined 
(IXiTj) . (jXsT]) into eq. ([02|) . we obtain 



(A.53) 



This agrees with the vacuum part of eq. (|A.58p which was independently determined with 
the method of sections IA.lHA.3l 



A. 5. Results for all independent sum-integrals 

We list here the results, after Fourier transform (eq. (|2.5p ) and cut (eq. ()2.6p ). for the "master" 
structures of eqs. (|3.18p . (|3.25p ~ (j3.28p . The bosonic and fermionic cases differ from each other 
simply by the exchange na -H" —np; therefore, we only list the bosonic expressions here. The 
results read: 
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167r3 
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g + a; w — a; — "^I/JJ 

Note that the integral J^dqqnB{q) appearing frequently could be carried out analytically 
(cf. eq. (I4.12|) ). but many others cannot so we choose to display all thermal contributions in 
a uniform unintegrated form. 



Appendix B. Details of Hard Thermal Loop resummation 

We present here the details of the Hard Thermal Loop resummation, outlined and motivated 
in sec. 13. 4[ Concretely, the computations below amount to a derivation of eq. (14. 3p . 

B.l. Framework 

At u -^T the unresummed result of eq. (14. 2p is logarithmically infrared divergent if we try to 
take the intercept determining the momentum diffusion coefficient k, defined in eq. p.7p (cf. 
eq. (I4.14P ). In order to obtain a result that is valid in this regime as well, the perturbative se- 
ries needs to be reorganized. The well-known framework for achieving this is that of the Hard 
Thermal Loop (HTL) effective theory |17j . In order to avoid double-counting, we first need 
to subtract the unresummed version of the HTL contribution from the QCD result, whereby 
we should obtain a result that is infrared finite but not physical, since the contribution of 
the soft modes, k ~ gT, is not included (cf. eq. (|2.9p ): subsequently the resummed HTL con- 
tribution is added, whereby the result remains infrared finite and now involves the physical 
contribution of the Debye scale as well. For both purposes, we start by recalling the HTL 
form of the gluon self-energy in D dimensions (at 0{g^) only the self-energy contribution is 
needed within the HTL theory, cf. section 3 of ref. [1]). 

Carrying out the Matsubara sums over g„ and in the self-energy of eq. (j3.19p : ex- 
panding the result to leading order in ikn/T, k/T; and restricting for the moment to spatial 
components, one finds 

U,,{K) =g\D-2)j^^ [2jVf n,{q) + {D - 2)N, n^iq)] + g'O , (B.l) 

where v = q/g. Introducing the independent projection operators 

P^AK) ^ 6,^5^, - Ml^ , pE^^K) ^ V - ^ - P^,AK) , (B.2) 
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and writing Iljj = P^^T + Pi^^E, the scalar structures Ht, can be projected out: 



Ue{K) = g\D-2 



[2Nin^{q) + {D-2)N^n^{q) k-v 



i/Cn + k • V 



Jlj,(K) = g2 f mnAq) + {D-2)N,n,{q) / k-v \ 



For future reference, we also introduce the notation 

1 r 



mi = g'{D-2) 



2N{n^{q) + {D ~ 2)N,n^{q) 



(B.3) 
(B.4) 



(B.5) 



which for D = 4 reproduces the known expression for the Debye mass parameter squared, 
given in eq. (j4.4p . 



B.2. Naive HTL analysis 

We first carry out an unresummed HTL computation, which needs to be subtracted from the 
naive QCD result in order to avoid double-counting and in order to obtain an infrared-finite 
result for the "matching coefficient" part of the final expression (cf. eq. (12. 9p ). This means 
that the self-energy is treated as an "insertion" , i.e. as if it were an interaction. We also note 
that since the HTL theory only represents the infrared physics of QCD and has not been 
derived with relative accuracy 0{g^), the gauge coupling appearing in it is to be treated as 
the renormalized coupling. Then the unresummed result contains a free contribution like in 
eq. dlS]), 
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6tt 
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as well as a next-to-leading order term which we first write in configuration space, 
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3 f^{K^f 



{D-2)klliT{K)+K^ IIe{K) 
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To arrive at this expression we started from the form in eq. ()3.20p and inserted = 
Hj- + P^H^;. Plugging in Hj-, H^; from eqs. ()B.3p . ()B.4p and partial fractioning the 
dependence on we get 
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It remains to Fourier transform (cf. eq. (j2.5p ) and take the discontinuity (cf. eq. (j2.6p ). For 
the structures at the end of the first row of eq. (jB.SP this goes Uke before (cf. eq. ()A.2ip ) and 
yields, for Z) — )• 4, 



K 



2 -^^^ A>7r 7~ 
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8^ 
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8^ 
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For the terms on the second row of eq. (|B.8p . on the other hand, the corresponding steps 
amount to 



ikn + A 



-k6{uj - A) 
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but the remaining integration over the spatial components of k requires regularization, so we 
have to proceed carefully. 

Making use of the integration measure in eq. ()A.17p . save for k, with the polar axis chosen 
in the direction of v = q/g and denoting z = k • v/A;, we are lead to consider the integral 
{d = D-l) 
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where a term with 5{uj + k) was omitted because of a restriction to > 0. In the term 
with 5{ijj — k) we can trivially integrate over k\ in the terms with 5{u] — kz) we can trivially 
integrate over z, noting that the constraints can get realized only for k > u. Renaming the 
integration variable subsequently as A: = coz, z > 1, and setting c? = 3 — 2e, the result can be 
expressed as 
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The first integral here can be rewritten as J^dy{y — l)"^"*^ with y = z^, and as an analytic 
function it vanishes. The second integral inside the square brackets evaluates to 1/e — ln4 + 
0(e). Expanding the prefactor as well and introducing the MS scale parameter as before, we 
obtain 



OJfl 



-2e / I 



Svr 



e 4a;2 



(B.13) 



Returning to eq. (jB.SP and inserting m| from eq. (|B.5P as well as the simple structures from 
eq. (|B.9p . we finally arrive at 
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B.3. Resummed intercept 

The remaining step is to carry out the resummed HTL computation (the last term in eq. (12. 9p ). 
As it happens, this problem was already addressed in ref. albeit for the difference 
p^{(jj)/uj — \m\^^^Q p^{(jj) /bj; the result is reproduced in eq. (|4.6p . In the following, it is 
then our task to consider the intercept Vrni^^Q p^{uj)/uj, which was not addressed in ref. [4]. 
In the following we define [pe]ir to be the spectral function computed within the HTL theory, 
but only up to linear order in a Taylor expansion in u around zero, i.e. corresponding to the 
intercept liiaii^^Q Pe{uj)/uj. It turns out that this limit is logarithmically ultraviolet divergent 
within the HTL theory, with the same 1/e as in eq. ()B.14p . 

The computation proceeds somewhat analogously to the unresummed one above, only now 
eqs. ()B.6p . (|B.7p . as well as higher order corrections are resummed. Writing the expression 
after the Fourier transform in eq. (|2.5p . and making use of the evenness of IIt,IIe in — ^ 
-kn,we get 



HTL, rcsum 3 



+ 



Taking now the discontinuity (eq. (j2.6p ) and noting that the first term of eq. (jB.lSP involves 
a high power of w, the contribution to the part linear in lo must come from the second 
term. Recalling eq. ()B.3p . we observe that + k'^ nicely cancels out. For small co we can 
furthermore write the structure inside H^; as 

^ '^ 1 + f( ] + iu}TT5{uj - k ■ v) ^ 1 + iujTr5{k ■ v) , (B.16) 



— + k • V — «0+ \— w + k-v^ 

where in the real part the principal value integral leads to a contribution vanishing faster for 
a; — )• than the unity. The unity leads directly to the m| of eq. (jB.SP and so we get 

-^(4)^^^1 ^ g^CpiP - 2)uj r 2N,n^{q) + {D-2)N,n^{q) f 7:kH{k • v) 

J IR, rcsum 

The k-integral here can again be carried out by choosing v as the polar axis and making use 
of the measure in eq. (|A.17P : 

d'^k fe2<5(k-v) 



vr 



(27r)'^ (fe2 + ml )2 

1 /■°° k'^^ 



(4^)^r(^bi) ^0 -(A;2 + m2)2 



-1 



d-3 



d^779 ^ dz {I - z') — 5{k z 



1 



2(47r)i-r(l - e) io (y + 1 
where y = k'^ /m\. Expanding in e we obtain 

Svr V e rn 



E 



^iR. = ^(T + lnf3-l) , (B.19) 
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and eq. ()B.17p in combination with eq. ()B.5P then yields 

r (4). J (fCpm^uHi^f I \ , . 

Pe^ = 5 s - + -1 • vB.20 

L J iR, losum 3 OTT \ e mg / 

Combining (as dictated by eq. (|2.9p ') the unresummed HTL result (eqs. (|B.6p . (|B.14p ) with 
the resummed IR limit (eq. (|B.20p ) we arrive at eq. (|4.3p . from which 1/e-divergences have 
duly cancelled. 

B.4. Resummed tj-dependence 

The remaining ingredient, the HTL-resummed contribution beyond the linear term in w, i.e. 
[Pe('^)]htl, rosum ~ ^e(^)]ij, ^^^^^ , was dctcrmiued in ref. g], and is reproduced in eq. (gS]) of 
the main text. The four terms correspond to the transverse cut, electric cut, transverse pole, 
and electric pole, respectively. The w-independent subtraction in the electric cut corresponds 
to eq. (fBTTll . 
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